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Abstract. Jet fragmentation functions provide insight into jet structure and 
are expected to be modified by the nuclear medium in A+A collisions with 
^ ' respect to p+p reference measurements. If jet reconstruction is unbiased then 

a softening of the fragmentation functions is expected and should be observed 
O ' in Au+Au collisions at RHIC. In these proceedings we present measurements 

2 , of fragmentation functions in p+p for charged particles for different jet find- 

ing algorithms; these measurements are understood and therefore can be used 
as a reference for comparison with Au+Au results. We report the effect of 
background and its fluctuations on jet reconstruction in Au+Au collisions, es- 
timated by using the jet algorithms on simulated Pythia jets embedded in real 
Au+Au events. Finally, measurements of fragmentation functions for jets re- 
, constructed in Au+Au events and their comparison to the p+p baseline are 

Tlj- ' presented and discussed. 
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H ' 1- Introduction 

, 

High-pT physics via di-hadron correlations has been a successful tool to probe the 
dense medium created in heavy-ion collisions at RHIC. Nevertheless, well-known 
geometrical biases [[I] provide limitations to such measurements. Full reconstruction 
of jets in heavy-ion collisions should allow direct access to the parton kinematics 
independent of the fragmentation process in an unbiased way. The first pt spectrum 
of jets reconstructed with a cone algorithm in p+p collisions at RHIC energies was 
measured by STAR with the 2003-2004 data, showing good agreement with NLO 
pQCD calculations over the measured range up to 50 GeV [El j^. Recent results 
from STAR [T, '5] show the feasibility of full jet reconstruction at RHIC in the high- 
multiplicity environment of Au+Au collisions. The underlying background and its 
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large fluctuations on an event-by-event basis make these measurements a challenge 
in the high- multiplicity environment at RHIC. Jet fragmentation functions provide 
insight into the structure of jets and are expected to be modified by the nuclear 
medium with respect to p-t-p reference measurements [ 6]. An expected softening 
of the fragmentation functions, i.e. a distortion of the so called "hump-backed" 
plateau [|7], should be observable in Au-|-Au with an unbiased jet population. 

2. Technical aspects in p+p and Au+Au analyses 
2.1. Event selection 

The STAR detectors employed for jet analyses are the Time Projection Cham- 
ber (TPC) that reconstructs charged particles and the Barrel Electromagnetic 
Calorimeter (BEMC) for the neutral energy. Corrections for double-counting of 
electrons and hadronic energy deposition in the BEMC are applied. The following 
trigger setups were used: 

(i) Jet Patch Trigger (JP) in p-l-p: year 6 run, with the transverse energy in a 
BEMC cluster Ary x Acf) — 1x1 Et >8 GeV (in addition to p+p minimum bias 
condition). The total luminosity is 8.7 pb^^. 

(ii) High Tower Trigger (HT) in p-fp: year 6 run, transverse energy in a single 
BEMC tower Arj x = 0.05 x 0.05 Et >5.4 GeV (in addition to p-fp minimum 
bias condition). The total luminosity is 11 pb^^. 

(iii) Minimum Bias Trigger (MB) in Au+Au: year 7 run, with coincidence condition 
between the two ZDCs. 

(iv) High Tower Trigger (HT) in Au+Au: year 7 run, transverse energy in a BEMC 
cluster A77 X A0 = 0.1 x 0.1 Et >7.8 GeV (in addition to MB condition). The total 
luminosity is 500 iJh~^, corresponding to a p+p equivalent of 19.6 pb^^. 



2.2. Jet-finding algorithms 

Jet-Finding algorithms cluster charged tracks and neutral towers into jets, which 
are experimentally defined as sprays of coUimated particles in the rj — cj) plane. Cone 
and recombination algorithms (see Ref. [[8] for more details) are the two main jet- 
finding tools. Cone algorithms start from those bins ("seeds") in the rj — (/) plane 
whose corresponding energy (from tracks and towers) is above a given threshold. 
For each seed, the algorithm searches for other particles in a cone of radius R 
{R = V(A(/>2 + A?72)). Cone jet-finders can be optimized by introducing iteration 
procedures and splitting/merging techniques [[5]. The infrared safety is an issue 
in jet algorithms and is partially solved by adding extra midpoint seeds. A proper 
solution is to use a seedless cone algorithm, like SISCone (Seedless Infrared Safe 
Cone) [[5], which is part of the Fast Jet [[TU] package. Recombination algorithms are 
seedless and cluster particles according to their spacial separation with a distance 
scale R and weight criteria. They are not generally bound to a circular structure. 
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Fig. 1. Uncorrected ^ distributions for four different jet algorithms for 30 < prjet < 
40 GeV with JP data, R=0.4 (left) and R=0.7 (right). 



In the "fcr" algorithm the weight is given by the particle pr and therefore soft 
particles will tend to cluster among themselves first. In contrast, in the "anti-Zcr" 
algorithm [ [11] the distance is weighted by the inverse px and hard particles are 
clustered together long before soft particles. Recombination algorithms are coUinear 
and infrared safe. 

3. Fragmentation functions in p+p 

The variable we use to describe the fragmentation function is ^ = lii(pT,jet/PT,h), 
where px.jet is the reconstructed transverse momentum of the jet and px.h is the 
transverse momentum of the charged hadron in the jet. As an alternative to z = 
PT,h/pT,jet, C is used to focus on the low-p^ part of the fragmentation function. We 
measured fragmentation functions from p+p JP data using four jet-finder algorithms 
(Seeded cone, "fcr" , "anti-/cT" , SISCone) for jet energies ranging from 10 GeV to 40 
GeV and above and for two choices of the jet radius, R=0.4 and 0.7 [[121 [ISI- It is 
known from Tevatron data that almost all the jet energy is contained in R=0.7 and 
more than 70%, increasing with the jet energy, lies within R<0.3. R=0.4 was used 
for the purpose of comparing p+p fragmentation functions with Au+Au results, 
where a small jet radius is necessary to suppress the huge background. Fig. [T] shows 
^ distributions for 30 < pTjet < 40 GeV with R=0.4 (left) and R=0.7 (right). 
Despite conceptual differences in the jet-finding algorithms studied, they agree over 
the range of jet energies explored. The observation that jet algorithms show similar 
performance with different jet radii is an indication that second-order effects like 
out-of-cone radiation are minor at RHIG energies, otherwise jet algorithms could 
behave differently especially in case of small jet radii. 
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4. Au+Au results 

4.I. Background effect in jet-finding in Au+Au 

As already mentioned, the background is the most crucial aspect in jet-finding in 
heavy- ion collisions. Possible ways to suppress the background can be to reduce 
the jet radius and require a minimum transverse momentum (p^"*) for a track or a 
BEMC tower to be included in the jet. On an event-by-event basis, the background 
contribution in the jet area is calculated as the mean pT in out-of-cone areas in the 
rj ~ (j) plane and is of the order of 45 GeV for R=0.4 and no Pj"*. In addition to 
its magnitude, that is subtracted, the background exhibits significant fluctuations 
in ?7 — regions on an event-by-event basis. Background fluctuations in these data 
can be estimated using a Gaussian background model, with the width decreasing 
from 6 GeV with no pf"* to 1.5 GeV with p™* = 2 GeV in areas of R=0.4 (Fig. El 
left). Background fluctuations in Au+Au, on top of the steeply falling jet spectrum, 
result in a shift towards higher energies of the jet px spectrum with respect to p-f p. 
This is shown in Fig. [5] (right), where spectra of Pythia jets reconstructed with a 
seeded cone algorithm (dashed lines) are compared to spectra obtained by running 
the same jet-finding algorithm on Pythia (version 8.1) jets embedded in real Au-|-Au 
events (solid lines). The effect of background fluctuations is dramatic with no pf"* 
(black lines) and is substantially reduced by applying a p™* — 2 GeV (red lines). 
This observation has important consequences in fragmentation function analyses, 
especially in view of comparing p-f p and Au-fAu measurements, where similar jet 
populations have to be selected in order to investigate possible modifications of 
the fragmentation functions. The most practical solution is to apply a p™* = 2 
GeV on both p+p and Au-|-Au events, that should result in a selection of jets with 
similar energies in p+p and Au+Au based on Pythia studies. The results of this 
analysis are reported in the next section. Alternatively, it is possible to avoid a p™* 
by unfolding the background fluctuations from the fragmentation functions and 
extract a distribution that is directly comparable to p+p. This work is in progress. 

4-2. Fragmentation functions in Au+Au 

The jet energy is determined in R=0.4 and with a p™* to suppress the background 
and its fluctuations. In order to include also the soft jet fragments, the fragmenta- 
tion functions are measured from charged hadrons in a radius of 0.7 around the jet 
axis without any minimum cut on transverse momentum of tracks. The measured 
fragmentation functions contain both jet and background components. The latter 
significantly contaminates the low momentum region of the fragmentation function 
and has to be subtracted. The background component of the fragmentation function 
is estimated on an event-by-event basis from charged particle spectra in the out-of- 
cone area. As shown in Fig. [3] (left), the fragmentation functions of mono-energetic 
30 GeV Pythia jets with and without embedding in Au+Au events are in agree- 
ment modulo small deviations due to the jet energy resolution: this suggests that the 
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Fig. 2. Left plot: width of background fluctuations as a function of jet radius 
for three choices of p™*. Right plot: reconstructed Pythia jet spectrum with and 
without embedding in 0-20% MB Au+Au events at \/s^^ = 200 GeV. 

background subtraction procedure provides a reliable measurement especially of the 
low-^ region of the fragmentation functions assuming Pythia fragmentation. The 
measured ^ distributions in p+p and Au+Au HT events are shown in Fig. [3] (right), 
for > 30 GeV. The corresponding jet momentum that selects similar jet 

populations in Au+Au as in p+p is pljl,f{AuAu) > 31 GeV for p™* — 2 GeV and 
was estimated with Pythia simulations (as described in the previous session). The 
background contribution from Au+Au jets with piy^^{AuAu) > 31 GeV (Fig. [3l 
right) starts to be dominant for ^ > 2.5 — 3. The ratio of the above fragmentation 
functions is shown in Fig. 2] for two choices of pj"* . The ratio close to 1 may sug- 
gest that the fragmentation functions in Au+Au HT events are not modified with 
respect to the reference p+p HT in the range 0.5 ^ ^ ^2.5. Further investigations 
are needed to understand the full implication of the result. 

5. Discussion of the results 

Results on jet spectra in Au+Au MB ([[S],[[T3]) suggest a binary scaling of jet 
production with respect to p+p within large systematic errors. With an unbiased 
jet population a softening of the measured fragmentation functions in the nuclear 
medium would be expected in order to explain the high-p^ hadron suppression 
observed at RHIC [ [151 Ull ■ On the contrary, no apparent modification of frag- 
mentation functions is observed in Au+Au HT events. Therefore, use of the HT 
trigger together with a p^"* seems to select a biased jet population. Two effects 
may occur: first, the HT trigger in Au+Au may select jets that are simflar to p+p 
because of "surface" effect and therefore are not modified by the medium. Second, 
the pj"* = 2 GeV (chosen because similar jet populations are selected in p+p and 
Au+Au assuming Pythia fragmentation) could largely affect the measurement if 
jets are softened in the medium and as a result the reconstructed energy would be 
underestimated . 
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Fig. 3. Left plot: ^ distributions of 30 GcV Pythia jets (open symbols), 30 GeV 
Pythia jets embedded in 0-20% MB Au+Au events with background subtraction 
(close symbols) and no background subtraction (dotted histogram). Right plot: 
uncorrected ^ distributions measured in p+p HT (open symbols) and 0-20% Au-f Au 
HT events at ^s^^ = 200 GeV for pl%t{PP) > 30 GeV and p™* = 2 GeV, with 
background subtraction (close symbols) and no background subtraction (dotted 
histogram). Only statistical errors. 




Fig. 4. Ratio of the ^ distributions measured in 0-20% Au+Au HT events at 
y/sj^ii^ — 200 GeV to p+p HT events at y/s = 200 GeV. The shaded area indicates 
the ^ region where the uncertainties due to background subtraction are large. 
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6. Conclusions and outlook 

The measurements reported in these proceedings show that fuh jet reconstruction 
is feasible at RHIC, although especially challenging in the heavy-ion environment. 
The jet fragmentation functions were measured in p+p with different jet-finding 
algorithms. The measurements provide a solid reference for Au-|-Au; they also sug- 
gest that out-of-cone radiation seems to be small at RHIC energies for cone radii 
R=0.4 and 0.7. The fragmentation functions in Au-|-Au for jet energies of 30 GeV 
and above were reported. Background fluctuations in Au-|-Au were reduced with 
a on particles. No modification was observed in the ^ distribution, indicating 
that a biased jet population was selected. A possible way to avoid the trigger bias 
is to study di-jets in the Au-|-Au HT events and measure fragmentation functions 
of recoil jets with respect to the trigger jets. In order to avoid biases due to the 
Pj"*, the analysis can be done with no p™* on the recoil jet, but in this case back- 
ground fluctuations will need to be unfolded to get fragmentation functions that 
can be directly compared to p-l-p. In addition, given the complexity of such mea- 
surements, data driven correction scheme and systematic error estimates (fake jets, 
detector effects, etc) are needed to have more control on our understanding of the 
jet fragmentation functions in Au+Au. These analyses are currently in progress. 
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